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Symmetry-protected topological (SPT) phases of matter have been the focus of many recent 
theoretical investigations, but controlled mechanisms for engineering them have so far been elusive. 
In this work, we demonstrate that by driving interacting spin systems periodically in time and 
tuning the available parameters, one can realize lattice models for bosonic SPT phases in the limit 
where the driving frequency is large. We provide concrete examples of this construction in one 
and two dimensions, and discuss signatures of these phases in stroboscopic measurements of local 
observables. 


I. INTRODUCTION 

Since the discovery of the quantum Hall effect 
(QHE) [1], topological phenomena in quantum rnany- 
body systems have dramatically changed our under¬ 
standing of phases of matter. In particular, the study of 
the fractional QHE brought about the notion of topo¬ 
logical order [2-4], which characterizes phases of mat¬ 
ter with emergent fractional excitations and topologi¬ 
cal ground-state degeneracy, which cannot be described 
within the standard Landau-Ginzburg framework. 

In recent years, the prediction and discovery of topo¬ 
logical band insulators [5, 6] has awakened a great deal 
of interest in gapped symmetry-protected, topological 
(SPT) phases of matter. These phases of matter lack 
fractionalized degrees of freedom, but display topolog¬ 
ical properties that manifest themselves in non-trivial 
boundary states that are protected by global symme¬ 
tries. While they do not display the long-range entan¬ 
glement of topologically-ordered systems, SPT phases of 
matter are characterized primarily by a nontrivial short- 
range entanglement structure in the low-energy states 

[7]- 

Following the classification of weakly-interacting 
fermionic SPT states [8-10], there has been a vast 
amount of recent effort to classify strongly-interacting 
SPT phases [7, 11-15] as well as to construct models 
supporting them [7, 16-25]. In light of this effort, it is 
highly desirable to identify controlled mechanisms ca¬ 
pable of bringing SPT states into realization. 

In this paper, we put forward a proposal to real¬ 
ize bosonic SPT phases as out-of-equilibrium states of 
quantum spin systems with periodically-driven multi¬ 
spin interactions. The systems we study are described 
by time-dependent Hamiltonians of the form 

H{t) =H 0 + 0(f) /(f) Hi nt , (1) 

where Hq is a local Hamiltonian describing a trivial 
paramagnet (i.e., one whose ground state is a trivial 
product state) and H lnt is a local interaction with a 
time-periodic coupling constant /(f) = /(f + T) with 


zero mean and a characteristic frequency w = 2n/T. 
0(f) is the Heaviside function denoting a protocol where 
the drive is switched on at t = 0. 

When Hi n t = 0, H(t) = Ho can be mapped from a 
trivial paramagnetic Hamiltonian to an SPT Hamilto¬ 
nian by a product of local unitary transformations that 
entangles the local degrees of freedom in a nontrivial 
way [7, 21]. Such transformations arise naturally in the 
study of many-body systems with periodically-driven in¬ 
teractions. In particular, we will show that, in the limit 
of large w, the time-periodic unitary transformation to 
the “rotating frame,” (we set h = 1) 

Unit) = e [ dt ' /(*')= e i9(t) , ( 2 ) 

generates the desired entanglement if H lnt is chosen ap¬ 
propriately. The transformation Unit) maps a state 
| ), whose time evolution is governed by the Hamil¬ 

tonian (1), into a state | ip R (t )) = U R (t) \ 4>{t)) whose 
time evolution is generated by 

Hn(t) = U n (t) H(t) t4(t) - i Unit) d t U^{t). (3) 

The stroboscopic evolution of the initial state in the ro¬ 
tating frame, \i() R {nT)) = e _1nT | iPr(0 )) (n G Z), 
is governed by the Floquet Hamiltonian Hf, which 
can be systematically determined via a Magnus expan¬ 
sion [26, 27]. (Note that Hf is also the generator of 
stroboscopic evolution in the “lab frame,” although we 
work with states in the rotating frame for convenience.) 
In the infinite-frequency limit, the Floquet Hamiltonian 
is nothing but the time-average of Hjt(f), 

H^ 0) = 2 j\tHn(t), (4) 

while the n-th order term in the Magnus expansion is 
of order l/w n . We will refer to Hp°^ as the stroboscopic 
Hamiltonian, because in the infinite-frequency limit, 
where only Hp°^ survives, the stroboscopic evolution and 
the true unitary evolution of the time-dependent system 
coincide. 
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FIG. 1. (Color online) Couplings in the leading term of the 
Magnus expansion (4) as functions of the scaled driving am¬ 
plitude A. White and gray regions correspond, respectively, 
to trivial and stroboscopic SPT phases. 

If the amplitude of the drive is small compared to the 
frequency, the stroboscopic Hamiltonian (4) simply re¬ 
duces to Hq. On the other hand, when the amplitude 
of the drive is chosen to scale with the frequency ui, the 
stroboscopic Hamiltonian can acquire a nontrivial form 
that is different from Hq [27]. In this work, we show 
that the stroboscopic Hamiltonian (4) describes micro¬ 
scopic models of SPT states with Z 2 x Z 2 [7, 20, 21] and 
Z 2 [7, 17] symmetries, respectively, for one- and two- 
dimensional driven systems. We refer to the phases gen¬ 
erated in this way as stroboscopic SPT (SSPT) phases. 
Remarkably, we find that, while the SSPT Hamilto¬ 
nian (4) is invariant under the global symmetry, the 
original time-dependent Hamiltonian (1) is not. Hence 
the global symmetry of the SSPT phase is found to 
be an emergent property of the high-frequency limit of 
%y- These results can be generalized to other symmetry 
classes. Finally, we also demonstrate that the dynamics 
of local observables at stroboscopic times can be used to 
probe the nontrivial edge states of SSPT systems with¬ 
out the need to prepare the system in the ground state 
of T~Lf- 

II. Z 2 x Z 2 SSPT PHASE IN ID 
A. The model 

We begin by studying an open ID chain with N sites 
described by the time-dependent Hamiltonian (1) with 

N N-l 

H m (t) =hJ2 a i+ 0 O) /(*) a i a i+l > ( 5 ) 

i—1 i—1 

and f(t) = Aw cos (cot + ft (A > 0). Note that 
the driving amplitude is taken to scale linearly with 
the frequency, so that A is dimensionless. The Pauli 
operators erf (a = x, y, z) satisfy the onsite algebra 


[erf, of] = 2i Sij e a t, c of and the anticommutation rela¬ 
tion {erf, of} = 2 S a b- Furthermore, notice that the 
Hamiltonian (5) has an onsite Z 2 spin flip symmetry 
generated by S = JX (=1 a i- Henceforth, we set the en¬ 
ergy scale h = 1, with the understanding that the limit 
ui —> oo corresponds to taking u>/i. 

Upon making the transformation to the rotating 
frame, we find that, for ip = 0, the stroboscopic Hamil¬ 
tonian (4) is given by 

Jft 2 xZ 2 = Jo(2A) (of + O'}}) 

( 6 ) 

+ E [a(A)uf-6(A)af_X^ + i] , 

i =2 

where a(A) = \ [1 + J 0 (4A)], 6(A) = 1 — a(A) and J 0 (ft 
is the Bessel function of the first kind. 

Observe that the Hamiltonian (6) possesses a global 
Z 2 x Z 2 symmetry generated by S even = J], eV en °? and 
<Sodd = Ui odd corresponding to independent spin 
flips on the even and odd sublattices. However, the 
time-dependent Hamiltonian (5) has a Z 2 symmetry, 
rather than a Z 2 x Z 2 symmetry - in other words, this 
enlarged symmetry group is an emergent property of the 
high-frequency limit ui —> oo, as it appears only upon 
taking the time average Eq. (4). 

We plot the couplings a(A), 6(A), and Jo(2A) in Fig. 1. 
By varying A, one can tune the couplings such that 
a(A) > 6(A) or vice versa. We will argue below that 
the values of A for which a(A) = 6(A) are critical points 
of the effective Hamiltonian that separate a trivial in¬ 
sulating phase from an SPT phase. 

To continue our analysis of this Z 2 x Z 2 -symmetric 
Hamiltonian, it is instructive to rewrite it in terms 
of Majorana operators ai = (U,<i of)erf and ft = 
icq erf, [28, 29] which are Hermitian and satisfy the 
usual fermionic algebra. In terms of these operators, 
the Hamiltonian (6) reads 


”Hz2xz 2 = — i Jo(2A)(«i ft + ajv ftv) 



“i—1 Pi -1 “i Pi “»+1 Pi+l a i+ 2 Pi+2 


FIG. 2. (Color online) Competing dimerization patterns in 
the Z 2 x Z 2 SPT chain. Dotted and solid lines account, 
respectively, for the dominant dimerization patterns in the 
trivial and SPT phases. When only the 3-spin term is present 
in the model, a dangling spin is localized on the edges. 
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This Hamiltonian contains two types of terms that can 
be thought of as projectors onto two distinct dimeriza¬ 
tion patterns that encode the entanglement structure of 
the ground-state wavefunction (see Fig. 2). The pattern 
encoded by the ct,;/3, terms involves Majorana dimers on 
each site. It is “trivial” in the sense that, for a finite 
chain, the pattern pairs all Majorana operators. On the 
other hand, the /3j_iQj_|_i terms encode dimerization be¬ 
tween next-neighbor Majoranas of opposite types. This 
pattern is “nontrivial” in the sense that it leaves two un¬ 
paired Majoranas at each end of a finite chain, yielding 
a fourfold ground-state degeneracy as a signature of the 
entanglement structure of the SPT phase. Equivalently, 
one can see this fourfold degeneracy from Eq. ( 6 ), as the 
operators af and afij commute with the Hamiltonian if 
only the three-spin interaction contributes. 

The transition between these two patterns and the as¬ 
sociated phases occurs at the point a( A) = b( A), where 
the bulk gap closes. To see this, we combine the Majo¬ 
rana operators into complex fermions c\ = (a* + i/3j)/2, 
in terms of which the Hamiltonian becomes 

N -1 

^Z 2 XZ 2 ^ ^ 1 T H.C.) 

i=2 ( 8 ) 

N-l v ' 

+ M C i C i + M<4 C 1 + C N C N ) 

i—2 

where we have dropped a constant term, and where we 
have defined t = 6 (A), p = 2a(A), and po = 2J 0 (2A). 
This model is nothing but two decoupled copies (one on 
the even sublattice and one on the odd) of the Kitaev 
model for a ID p- wave superconductor [28]. The critical 
point for this model is well-known, and occurs for p = 
2 1. However, this is equivalent to the condition a(A) = 
6 (A). It is important to note that while the chemical 
potential po at the ends of the chain is not equal to 
the bulk value /i, the location of the transition is not 
affected for a sufficiently long chain, as we have verified 
by exact diagonalization. 

The preceding discussion illustrates that the strobo¬ 
scopic Hamiltonian ( 6 ) is Z 2 x Z 2 -symmetric and con¬ 
tains one free parameter, A, that tunes the system across 
the transition between the Z 2 x Z 2 SSPT phase and 
the trivial paramagnetic phase. However, it is evi¬ 
dent from Fig. 1 that the coupling a(A) 7 ^ 0 for any 
A. Naively, then, it seems that one cannot access the 
“ideal” scenario where the operators af N commute with 
the Hamiltonian. Nevertheless, this is not the case, as 
the local field Jo( 2 A) at the ends of the chain vanishes 
identically if 2A is equal to a zero of the Bessel func¬ 
tion J 0 (see Fig. 1). In this case, the operators of and 
a z N identically commute with the Hamiltonian ( 6 ), and 
the system has an exact fourfold ground-state degen¬ 
eracy despite the presence of a transverse field in the 
bulk. Deviations from these special values of A split 


this degeneracy by an amount that decreases exponen¬ 
tially with system size, and the system remains in the 
SSPT phase so long as the bulk gap remains open. 

B. Signatures in stroboscopic dynamics 

1. Infinite-frequency limit 

So far, we have demonstrated that, in the limit w —7 
00 , the stroboscopic evolution of the periodically-driven 
spin chain of Eq. (5) is generated by the effective Hamil¬ 
tonian ( 6 ), for an appropriately-chosen driving protocol. 
However, it remains to be shown that the stroboscopic 
evolution generated by Eq. (6) yields telltale signatures 
of SPT physics in local measurements. To address this 
point, we consider the local spin expectation value 

(of(t)) = (Wole+^^of e - i ^xz 2 t |^ 0 ), ( 9 ) 

where Tq ) is some initial state. This quantity coin¬ 
cides with the true time evolution of the operator of in 
the limit w —> 00 , where the period T is infinitesimally 
small and t = nT (n £ Z) is approximately a continuous 
variable. For simplicity, we choose the initial state | 'Po ) 
to be a product of eigenstates of af on each site j, so 
that (a f(t)) is invariant under Un(t) ( i.e the unitary 
transformation to the rotating frame). 

The observable defined in Eq. (9) provides a clear 
signature of the edge states in the SSPT phase, even 
though the product state | \Po) is a highly out-of- 
equilibrium state with respect to the stroboscopic 
Hamiltonian (c.f. Ref. [30]). When Jo(2A) = 0, (erf(t)) 
and (crfj(t)) are independent of time, since af and a z N 
commute with the effective Hamiltonian "Hz 2 x z 2 - For 
i 7 ^ 1 or N, however, (of (t)) evolves quasi-periodically 
in time, with oscillations occurring on a timescale r c 
on the order of the inverse bulk energy gap of H^x z 2 
[see Fig. 3(a)]. If A is tuned slightly away from one 
of these special values but remains within the phase 
boundary, which amounts to adding a small transverse 
field (af + afj), then the end spins precess in the y-z 
plane on timescales much longer than r c , so that the 
bulk and boundary behavior can be distinguished. Be¬ 
cause the end spins were completely unconstrained be¬ 
fore the addition of the transverse field, this preces¬ 
sion, which is characteristic of a free spin placed in a 
magnetic field perpendicular to the quantization axis, 
manifests itself in oscillations of (af N (t)) about zero 
[see Fig. 3(b)]. In this way, (af N (t)) can be used to 
distinguish the free edge spins characteristic of SPT 
physics from spins that are “frozen” due to the pres¬ 
ence of a symmetry-breaking field. For example, if a 
small longitudinal field h z J2iLi is added to Eq. (5), 
thereby breaking the Z 2 x Z 2 symmetry of the strobo¬ 
scopic Hamiltonian, then one finds that the end spins 
no longer wrap the unit circle in the y-z plane as they 
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FIG. 3. (Color online) Stroboscopic evolution of (of(t)) for im 1,2, and 5 for an eight-site chain, where the initial state is 
chosen to be the product state | \to ) = | tt • • • t) ■ Panels (a)-(c) depict stroboscopic evolution at u> = oo [i.e. as defined in 
Eq. (9)], while panels (d)-(f) depict the exact stroboscopic evolution [i.e. as defined in Eq. (10a)] for uj = 100 h. In (a), 2A 
equals the second zero of the Bessel function Jo(x) in Fig. 1, and the longitudinal symmetry-breaking field h z = 0. In (b), 
A = 2.6 and h z = 0, while in (c), A = 2.6 and h z = 0.01 in units of the bare transverse field h. The plots in panels (d), (e), 
and (f) use the same parameters as the ones in panels (a), (b), and (c), respectively. 


precess when the transverse field is added at the ends of 
the chain. Consequently, (v* N (t)) no longer oscillates 
around zero, but around a nonzero value whose sign 
matches the sign of the longitudinal field [see Fig. 3(c)]. 
Thus, in the infinite-frequency limit, stroboscopic mea¬ 
surements of (erf (t)) are an effective dynamical probe of 
SPT physics at the boundary of the driven system. 

2. Finite-frequency corrections 

The preceding discussion is an accurate description 
of the driven system in the limit where the driving fre¬ 
quency u> is infinite. However, at any finite frequency, 
there are corrections to this behavior that become im¬ 
portant in the infinite-time limit, where the “error” 


due to these corrections can accumulate without bound. 
We now characterize the nature of these corrections 
(c./. Ref. [27]), and present arguments and numerical 
results showing that there is a window of time after the 
drive is switched on during which the signatures of the 
infinite-frequency SSPT phase can be observed at large 
but finite driving frequencies. 

As the driving frequency oj is decreased, two effects 
occur that lead to deviations from the infinite-frequency 
case discussed in the previous section. First, the stro¬ 
boscopic evolution of observables, as in Eq. (9), no 
longer coincides with the true time evolution of the sys¬ 
tem. In particular, expectation values of observables 
become dressed by intra-period effects that become sig¬ 
nificant if the system is not observed at stroboscopic 
times t n = nT for n £ Z [27]. However, the expectation 
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values of observables at stroboscopic times are still pre¬ 
dicted by the unitary evolution generated by the Flo- 
quet Hamiltonian T~Lp. For instance, the stroboscopic 
evolution of erf, which is given at infinite driving fre¬ 
quency by Eq. (9), becomes 

<*?(*„)> = (*o| e +iH ^ a* e- iH ^\* 0 ), ( 10 a) 

where 

n F = n { p ] + n ( F 1) + ... (iob) 

contains all finite-frequency corrections 'Hp'* ~ h ( h/uj) k 
to the infinite-frequency Floquet Hamiltonian 'Hp\ 

This brings us to the second effect, namely the fact 
that Tip acquires finite-frequency corrections appear¬ 
ing at orders 1/w and higher in the Magnus expan¬ 
sion. These corrections generically break the symme¬ 
try that protects the SSPT phase. Indeed, in Appendix 
A, we present the leading finite-frequency correction to 
the infinite-frequency Hamiltonian %z 2X z 2 that break 
the emergent Z 2 x Z 2 symmetry down to the Z 2 sym¬ 
metry of the original time-dependent Hamiltonian (5). 
We will now argue that these symmetry-breaking cor¬ 
rections are unimportant for the detection of the stro¬ 
boscopic signatures of the SSPT phase discussed in the 
previous section, provided that the driving frequency w 
is sufficiently large compared to the bare transverse field 
energy scale h. 

To see this, observe that at some stroboscopic time 
t n = nT, the fc-th order Magnus correction is only im¬ 
portant if 'Up' 1 t n is comparable in size (modulo 27 t) to a 

number of order one. Since Hy k> ~ h (h/uj) k , this means 
that one must wait a (stroboscopic) time 



in order for symmetry-breaking effects to begin to mani¬ 
fest themselves. Note that, in the limit w —> 00 , the time 
f* = ?z* T —>• 00 as well, so that the infinite-frequency 
limit manifests the enlarged Z 2 x Z 2 symmetry, as ex¬ 
pected. 

To support this argument in the context of the 
dynamical signatures of the SSPT phase discussed 
in the previous section, we have supplemented the 
infinite-frequency stroboscopic evolution [c.f. Eq. (9)] 
displayed in Fig. 3(a)-(c) with finite-frequency calcu¬ 
lations [c.f. Eq. (10a)] over a range of frequencies. 
The finite-frequency stroboscopic calculations were per¬ 
formed using the exact evolution operator over a pe¬ 
riod, determined by direct numerical integration of the 
Schrodinger equation. We find that a driving frequency 
u! ~ 100 h is sufficiently large to extract the dynami¬ 
cal information necessary in order to infer the existence 
of symmetry-protected edge states in the manner out¬ 
lined in the previous section [see Fig. 3(d)— (f)]. Differ¬ 
ences between these finite-frequency results and their 



FIG. 4. (Color online) Finite-frequency corrections to stro¬ 
boscopic evolution of (of(t)). The system size and initial 
state are the same as in Fig. 3. Panel (a) shows the evo¬ 
lution depicted in Fig. 3(e) over a longer time, so that the 
envelope timescale discussed in the main text is more appar¬ 
ent. Panel (b) uses the same parameters as Fig. 3(a) and 
(d), but at frequency u> = 10 h. 


infinite-frequency counterparts only begin to manifest 
themselves on timescales of order f* = 10 4 T, where, 
according to the scaling argument in the previous para¬ 
graph, the k = 1 correction begins to become important. 

The main difference that arises at large but finite fre¬ 
quencies is the appearance of an envelope timescale, 
much longer than t*, on which the end spins oscillate 
(see Fig. 4). For example, in Fig. 4(a), the dynamics 
shown in Fig. 3(e) is plotted out to a time t = 10f*, 
where the appearance of the envelope timescale is very 
clear. At the “magic” values of A where the end spins are 
frozen in the infinite-frequency limit [c.f. Fig. 3(a)], the 
envelope timescale also manifests itself in oscillations 
of the previously-frozen end spins due to the appear¬ 
ance of symmetry-breaking finite-frequency corrections 
at long times. This is already evident in Fig. 3(d), where 
oj = 100 h, and is exaggerated in Fig. 4(b), where the 
frequency has been lowered to 10 h. 

To summarize this discussion, we have argued in this 
section that the dynamical signatures of the SSPT phase 
discussed in the previous section are immune to finite- 
frequency corrections for a window of time whose size 
increases monotonically with the driving frequency and 
approaches infinity in the infinite-frequency limit. At 
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high but still reasonable frequencies w ~ 100 h , this time 
window is sufficiently large to extract these dynamical 
signatures before the corrections begin to take over. 


III. RATIONALE, GENERALIZATIONS, AND 
POSSIBLE EXPERIMENTS 

The rationale behind the ability to engineer the SSPT 
Hamiltonian ( 6 ) can be stated as follows. First, re¬ 
call that, when the driving vanishes, the mapping from 
the trivial phase to the SPT one can be achieved via 
a product of local unitary transformations [7]. In the 
Z 2 x Z 2 case, the generator of this transformation is 
proportional to the Ising interaction in Eq. (5) [21], On 
the other hand, in the driven system, the unitary trans¬ 
formation Uft(t) to the rotating frame is also generated 
by the Ising interaction. Consequently, at infinite fre¬ 
quency, we found parameter regimes in which this trans¬ 
formation effectively mapped a trivial paramagnet to an 
SPT one, and gave rise to an emergent Z 2 x Z 2 symme¬ 
try that is not shared by the time-dependent Hamilto¬ 
nian (5). 

The above discussion suggests a principle for obtain¬ 
ing an SSPT phase in a periodically-driven system: the 
interaction term in Eq. (1) should be chosen to be the 
generator of the unitary transformation connecting a 
trivial to an SPT system. In order to demonstrate that 
this stroboscopic approach to SPT phases applies be¬ 
yond the ID case discussed above, we consider a 2D 
system on the triangular lattice with a driven three-spin 
interaction, 

H 2 u(t) = ~ h Yl + 0 (*) /(*) zL °i a J ’ (12) 

3 (ijk) 

where the summation in the second term runs over all 
the triangles of the lattice and we assume the same form 
for f(t) as in Eq. (5). Interestingly, for A « 0.51 and 
ip « ±0.277t, we find that (see Appendix B) 

Hp' 1 « h e ff ^2 a j e * ? ■ ( 13 ) 

j 

Remarkably, the Hamiltonian (13), which involves up 
to seven-spin interactions, is the exactly-solvable model 
of a Z 2 SPT paramagnet studied by Levin and Gu in 
Ref. [17]. The model (13) can be obtained from a trivial 
paramagnet by a product of local unitary transforma¬ 
tions that each depend on three a~ spins (see Appendix 
B), which then justifies the need for a three-spin inter¬ 
action in ( 12 ). 

We close by commenting on possible experimental 
realizations of SSPT phases. Recent developments in 


quantum simulation with trapped ions [31-34] and su¬ 
perconducting quantum circuits [35, 36] have shown 
that it is possible to engineer tunable multi-spin interac¬ 
tions and transverse fields in a laboratory setting. These 
developments suggest the possibility that SSPT phases 
like the ones discussed in this paper could be realized in 
an experiment, if the appropriate sinusoidal drive can 
be implemented. The superconducting quantum circuit 
architecture described in Ref. [36] appears particularly 
well-suited to these purposes, as it was demonstrated 
in that work that the couplings between the supercon¬ 
ducting qubits in that system can be tuned dynamically. 
The periodic modulation of the interaction strength re¬ 
quired by our proposal is already feasible in that setup, 
making it an ideal candidate for a possible experimental 
realization. A thorough assessment of the suitability of 
this proposal for the experimental platforms mentioned 
above is necessary, but beyond the scope of this work. 

To summarize, we have shown in this work that, 
by adding appropriately chosen periodically-driven mul¬ 
tispin interactions to a trivial paramagnetic Hamilto¬ 
nian, it is possible to realize SPT phases in the high- 
frequency limit. We further illustrated via a ID ex¬ 
ample that the SPT phase can be probed with stro¬ 
boscopic measurements of local observables. We also 
characterized numerically (and analytically, in the Ap¬ 
pendices) the finite-frequency corrections to the pure 
infinite-frequency SSPT Hamiltonian, and found that 
driving frequencies of order a hundred times the charac¬ 
teristic bare energy scale of the problem are sufficient to 
observe signatures of the phase. Finally, we illustrated 
how this construction can be extended to higher dimen¬ 
sions and different symmetry classes, such as the above 
example of the Z 2 SPT phase in 2D. In future work, 
it would be interesting to further explore the possibil¬ 
ity of generating nontrivial patterns of entanglement by 
driving. The work presented here can also be used as 
a springboard to future progress in the development of 
out-of-equilibrium and non-eigenstate probes of topo¬ 
logical and symmetry-protected topological order. In 
particular, it would be interesting to determine, along 
the lines of Ref. [30], how the presence of disorder can 
protect SPT order in the novel context, explored here, 
of periodically driven systems. 
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Appendix A: One-Dimensional SSPT Hamiltonian 
1. Stroboscopic Hamiltonian 

We derive an effective Hamiltonian that encapsulates the stroboscopic dynamics generated by 

N N-l 


Hin(t) = hJ2 fJ i + ©(*) fit) Y a i a i +1 ’ 


(Ala) 


with 


f(t) = Aw cos (cut + p) , A > 0 . 


(Alb) 


To do this, we employ the time-dependent unitary transformation 

,t JV-i 


JV-l 


U R {t) = exp i / dt' &{t) f(t) Y a ! a t 


i +1 


2=1 


= exp 


Y a i aZ 


i +1 


2=1 


where g(t) = A [sin (u; t + y?) — siny?], which transforms the Hamiltonian to the rotating frame as follows: 


N 


H R (t) = U R {t)H 1D (t)Ul(t)-iU R {t)dtU' R (t) = hU R (t ) £>? ul(t). 


Explicitly we find 
N -1 


#rW = Y { cos 2 ( 2 <?(t)) erf — sin 2 ( 2 g(t)) of_ 1 0 'f 0 f +1 — cos( 2 g(i)) sin( 2 <?(t)) (of<rf +1 + cr|_ 1 (r| / )| 
2=2 

+ cos(2g(t)) (a? + a x N ) - sin( 2 5 (t)) (erf erf + ct^) . 

The time average of Eq. (A4) yields 


N-l 


^F 0) = E { a ( A ’ - 6 ( A ’ ¥>) - c(A, <p) ( crfcrf +1 + cr.f_ lC rf ) j 

2=2 

+ rf(A, </?) (erf + cr^r) - e(A, <£>) (erf <T^ + cr^.! <J V N ) , 
where the coefficients a(A, p), ..., e(A, </?) are given by 


1 /’ 27r 

i(A, p) = 1 - b(A, ¥?) = —/ dr cos 2 [2AG v (r)j , 
27T Jo 

1 f 27r 1 

c(A, = 2 n J Q dT 2 Sin [ 4AG >( T )] ’ 


(A2) 


(A3) 


(A4) 


(A5) 


(A 6 a) 


(A 6 b) 


1 l’ 2 ”’ 

d(X,p) = — / dr cos [2AG v (r)] 

27T Jo 


(A 6 c) 


1 f 27 " 

e(A, p) = — / dr sin [2AG v (r)j , (A 6 d) 

27t Jo 

where G v (r) = sin (t + p) — siny>. Observe that for the choice p = 0, the two-body terms that break the Z 2 x Z 2 
symmetry vanish and we recover the stroboscopic Hamiltonian in Eq. ( 6 ) of the main text. 
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FIG. 5. (Color online) Couplings in the first-order correction to the Magnus expansion, c.f. Eq. (A8), as functions of the 
scaled driving amplitude A. White and gray regions are used to distinguish the trivial and stroboscopic SPT phases, as in 
Fig. 1, and the Bessel function Jo(2A) is plotted for reference. 


2. Leading Finite-Frequency Correction to the SSPT Hamiltonian 


We now present the order-1/w Magnus correction to the SSPT Hamiltonian, namely 


U/(l) _ 

rip — 


2 T 


dti / dt 2 [H R (t 1 ),H R (t 2 )], 


(A7) 


where H R (t) is given by Eq. (3). We will work exclusively with an infinite chain in this section, as our aim is only 
to show that the bulk Z 2 x Z 2 symmetry is broken by this correction. After calculating the necessary commutators, 
we find that 


Wp =h 1 Y J (c 


'i +1 


- oM +1 ) +h 2 J 2 K-P 


T i+l a i+2 


where the coefficients are given by 


K 

h 2 



T 1 

dr 2 cos(2Asinri) cos(2Asinr 2 ) sin[2A(sinr 2 
dr 2 sin(2A sin Ti) sin(2Asinr 2 ) sin[2A(sinr 2 — 


- sinn)] 
sin n)]. 


(A8a) 


(A8b) 

(A8c) 


These coefficients are plotted as functions of A in Fig. 5. Observe that each term above breaks the Z 2 x Z 2 symmetry. 
The Z 2 x Z 2 symmetry of the zeroth-order Hamiltonian ( 6 ) is therefore an emergent symmetry that appears only 
at high frequencies. 
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FIG. 6 . Triangular lattice where the model Eq. (B2) and the driven Hamiltonian Eq. (B 6 a) are defined. The sites nearest 
neighbors to site j are labeled 1 through 6 . 


Appendix B: Two-Dimensional SSPT Hamiltonian 
1. Exactly Solvable SPT Model 

In this section, we review the 2D Z 2 SPT model introduced by Levin and Gu in Ref. [17]. We start with the 
trivial paramagnetic Hamiltonian on the triangular lattice (see Fig. 6 ), 




(Bl) 


and the 2D Z 2 SPT Hamiltonian [17] 

H2D,spt = — ^2 = ^2 a j e ' 4 ^ a ‘ > 


(B2) 


■3 0 


where the sum over £, £' in Eq. (B2) extends over pairs of nearest neighbor spins around the spin at site j as depicted 
in Fig. 6 . The Hamiltonian Eq. (B2) is invariant under spin flips generated by Sz 2 = J][ 

The Hamiltonians Eq. (Bl) and (B2) are related by the unitary transformation 

g _1 5 ( 1-<T J ) _1 f ( 1_<T I cr*,) ^ 

j 

that implements 


Oj = WtrJW 


_1 = — uj e ^ ( 1_ °7 <T I') 


(B3) 


(B4) 


Now for every site j we expand the exponent: 

g 1 4 52<«'>o CT «') = e 1 ^ [cos ( 7 r/ 4 ) — i sin ( 7 t/ 4 ) a\ erf,] 

_ ^ / ,^Z „Z „Z Z Z I ^.Z „Z Z Z I ^.Z Z Z Z I Z ,Z ^_Z Z 1 Z Z Z Z I ^.Z ,Z Z „Z I ^.Z ^.Z Z ,Z (B5) 

— ^ °2 a 3 a 4 a 5 a 6 + a l a 2 a 3 a 5 + a l a 3 a 4 a 5 + 17 1 a 2 a 4 a 6 + a 2 a 3 a 4 a 6 + <J 1 °3 ^5 °6 + CT 2 °4 °5 CT 6 

I —2: —Z | _21 _2; i -Z ^.Z (^Z ^.Z I _Z —Z i ^2 ^2 , _2 _2 i ^2 _2 I _2 2 \ 1 

+ (T 1 cr 4 + <T 2 a 5 + a 3 a 6 ~ \ a l a 2 + a 2 a 3 + a 3 °4 + °4 CT 5 + CT 5 °6 + °6 ° 1 j j ) 

where erf, ...,<rf denote the six spin operators around the site j, as in Fig. 6. 
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2. Driven Three-Spin Interaction 

Motivated by the unitary transformation Eq. (B3), we are led to consider a time dependent three-spin interaction 


H 2D (t) = -h Y erf + 0 (f) /(f) J2 °i a j a k > 

j ( ijk ) 

where the summation in the second term runs over every triangle of the lattice and 

/(f) = Aw cos(wf + p ), A > 0. 


(B 6 a) 


(B 6 b) 


The unitary transformation to the rotating frame 


Z7 R (t)=exp ig(t) Y ViVjVk 

{ijk) 


where g{t ) = A [sin(wf + ip) — sin yp] , yields the rotating-frame Hamiltonian 


(B7) 


H R (t) = U R (t) -h Y Ul(t). 


(B 8 ) 


The relevant object to compute is then 


Ui i(f) aj t 4 (t) = af exp 


Explicitly, we have 


-12 g(t)(Tj Y a e a 

(M'j) 

= erf n [ cos (2g(f)) - i Sin(2g(f)) erf erf erf,] 
(Wj) 

= a j A' (*) ■ 


-4j(f) = -4](f) + M“(f) 


where 


(f) = i a 2 [ft(f) (af af + ag erf + er: 


3 u 6 


1 _Z _Z _Z _Z I _Z _Z _Z _Z I _z _z _z _z 

+ a l a 2 a 3 a 5 + a l a 3 a 4 a 5 + °1 a 2 a 4 a 6 

1 _2 _Z „Z ^.Z I _Z ^Z ^.Z ^Z 1 _z ^z ^.z ^.z 

+ er 2 cr 3 cr 4 a 6 + a 1 a 3 a 5 a 6 + a 2 a 4 a 5 a 6 


a l a 2 a 3 a 4 a 5 a 6 


Q (j-\ („Z „Z I ^Z Z I _Z _Z I _Z _Z I _Z _Z I _z ^.z \ 

— P 2 (f) (er 4 er 2 + a 2 a 3 + a 3 cr 4 + cr 4 a 5 + cr 5 cr 6 + er 6 cr 4 ) 


(B9) 


(BlOa) 


(BlOb) 


A l /(t) = fait) + /3 4 (f) (af erf erf ag + af ai af af + af af af af 
+ erf af af ag + af ag af ag + af af ag ag 

1 —2 _Z _Z _Z I _Z ^.Z _Z _Z I _z _z _z _z 

+ <T 2 er 3 a 4 a 5 + a 2 a 3 a 5 a 6 + a 3 a 4 a 5 a 6 

I ^-2 _2 I _Z _Z I _Z _Z I _Z _Z I _Z _Z I _Z _z\ 

+ a i a 3 + <r 4 a 5 + a 2 a 4 + a 2 a 6 + a 3 a 5 + a 4 a 6 J , 


(BlOc) 
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where 

0i(t) = 2 c 3 (t) s 3 (t ), 

P 2 (t) = c(t) s 5 (t) + c 5 (t) s(t ), 
As(i) = c 6 (t) - s 6 (t ), 

/3 4 (i) = c 2 (t) s 4 (f) - c 4 {t) s 2 (t ), 


(BlOd) 


and we use the shorthand notation c(t) = cos( 2 g(t)) and s(t) = sin( 2 g(t)). 

The Floquet Hamiltonian at infinite frequency, obtained from the time average of the Hamiltonian Eq. (B 8 ), 


?d 0) = -h 


E 

3 


if l dt - 4 >«) 


fun) 


upon using Eq. (BIO), depends on the following parameters 


/?i (A, v?) 


1 f 2v 

— / dr 2 cos 3 2AG v (r) 

2 7T ,/n 


sin 


2A G v (t) 


(B12a) 


@2{X,ip) = — 

Z 7r 




cos 


2 A G*(t) 


sin 


2AG^(r) 


cos 


2 A G*(r) 


sin 


2AG„(t)] } 


(B12b) 


/?3(A, v?) 


1 

2 7T 



2AG ¥ ,(t) 


-sin 6 



(B12c) 


= — 
Z 7 r 




cos 


2AG< /? (r) 


sm 


2AG^(r) 


— cos 


2AG^(r) 


sm 


2AG v (t)] } , 


(B12d) 


where G v (t) = sin (r + ip) — sin 93 . Whenever 


ft(A*, V *) = ft(A‘,/)^V0, 


(B13a) 


/3 3 (A*,^) = /3 4 (A*,^) = 0, 
the Hamiltonian Eq. (Bll) acquires the form 

7^ =4 P* h ^2 e 1 * ^<“'>0 ( 1_cr ? <7 I') 

j 

which is the same model Eq. (B2) shown in Ref. [17] to describe the 2D SPT paramagnet with Z 2 symmetry. We 
have found numerically that condition Eq. (B13) is satisfied, for example, for A* ss 0.51 and ip « ±0.27 7 r. It is 
fundamental to stress that even though the driven Hamiltonian Eq. (B 6 a) does not have the Z 2 symmetry, this 
symmetry emerges in the u> —> 00 limit. 


(B13b) 

(B14) 
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